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ANALYSIS  OF  THE  IIP  ICEBERG  DRIFT  MODEL 


ABSTRACT 

The  International  Ice  Patrol  has  developed  a  model  of  iceberg  drift  that 
plays  a  critical  role  in  the  overall  modeling  of  iceberg  locations  and 
determining  the  Limits  of  All  Known  Ice.  The  drift  model,  along  with  an 
iceberg  deterioration  model,  is  used  to  determine  the  location  and  state  of 
any  reported  iceberg  until  it  is  resighted.  The  drift  model  depends  on  the 
reported  size  of  the  iceberg  and  a  number  of  environmental  parameters. 
This  report  examines  the  structure  of  the  drift  model,  reviews  empirical 
evaluations  of  the  model,  and  conducts  an  analytical  and  experimental 
sensitivity  analysis  of  the  model  outputs  with  respect  to  the  model 
parameters.  The  result  of  this  analysis  suggests  that  the  local  wind  driven 
current  portion  of  the  deterioration  model  should  be  reexamined.  The 
analyses  indicate  that  the  current  speed  magnitude  is  somewhat  sensitive  to 
errors  in  the  wind  speed  that  is  provided  as  an  input.  Similar  analyses  with 
respect  to  errors  in  wind  direction  indicate  that  errors  in  wind  direction 
may  have  a  significant  impact  on  estimated  current  direction.  A 
comparison  with  a  SAR  current  model  yielded  significantly  different 
results.  Before  any  current  model  changes  are  made  based  on  the  analysis 
of  the  SAR  current  model,  the  differences  between  the  two  models  must  be 
resolved.  Other  aspects  of  the  drift  model  appear  to  provide  a  reasonable 
representation  of  the  actual  drift  process.  The  analytical  evaluation  of  the 
iceberg  drift  model  reveals  that  there  is  little  need  for  improved  estimates 
of  the  environmental  parameters  for  direct  use  in  the  drift  model.  The 
analysis  illustrated  the  importance  of  correct  classification.  It  is  suggested 
that  the  current  policy  of  classifying  unknown  icebergs  as  non-tabular 
medium  icebergs  be  reexamined.  Finally,  the  sensitivity  of  positional 
accuracy  was  clearly  illustrated  for  the  drift  model.  It  is  important  to  be 
able  to  improve  the  initial  positioning  accuracy  to  ensure  that  the 
probability  of  using  the  correct  geostrophic  current  is  maximized. 


INTRODUCTION 


Objective. 

The  IIP  uses  an  Iceberg  Drift  Model  to  emulate  the  approximate  drift  behavior  of 
icebergs  in  the  presence  of  known  environmental  conditions.  The  results  of  this  model 
along  with  those  from  the  iceberg  deterioration  model  result  in  a  description  of  the 
location  and  melt  sate  of  icebergs  that  is  used  to  determine  the  Limits  of  All  Known  Ice 
(LAKI)  The  selected  modeling  alternatives  for  Phase  II  of  the  Cost  and  Operational 

-  nd  /  or 
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Analysis  of  the  IIP  Iceberg  Drift  Model 


Effectiveness  Analysis  included  conducting  a  detailed  sensitivity  analysis  of  the 
deterioration  and  drift  models  and  developing  an  approach  to  characterize  the  risk  posture 
for  the  IIP  (Armacost,  1994).  The  detailed  sensitivity  analysis  of  the  deterioration  model 
and  the  development  of  a  risk  model  are  addressed  in  other  reports.  The  purpose  of  this 
report  is  to  review  the  structure  of  the  existing  model  and  to  exanune  the  sensitivity  of  its 
output  with  respect  to  changes  in  input  parameters.  This  analysis  identifies  which  input 
parameters  require  the  most  attention  with  respect  to  accuracy  of  their  estimates  and 
identifies  areas  where  potential  model  enhancements  are  appropriate. 

Background. 

The  iceberg  drift  model  used  by  IIP  was  completed  and  tested  in  1980  (Mountain, 
1980).  Its  format  and  use  remains  essentially  the  same  to  date.  The  fundamental  drift 
model  balance  is  between  iceberg  acceleration,  air  and  water  drag,  the  Coriolis 
acceleration  and  a  sea  surface  slope  term  characterized  by  the  mean  ocean  currents.  The 
resulting  differential  equations  are  solved  using  a  fourth-order  Runge-Kutta  analysis.  Key 
input  data  include;  iceberg  location;  iceberg  size  and  shape;  local  wind;  geostrophic 
(mean)  current;  and  local  currents  from  drift  buoys. 

Sensed  or  sighted  icebergs  are  placed  into  one  of  four  size  categories  (growler, 
small,  medium,  and  large  with  no  upper  limit)  which  also  automatically  sets  the  mass  and 
cross-sectional  areas  to  the  assumed  characteristic  values  for  the  designated  size  category. 
One  of  two  specific  shape  classifications,  tabular  or  non-tabular,  is  also  made  when  a 
visual  sighting  occurs.  The  model  divides  the  subsurface  shape  of  the  iceberg  into  up  to 
four  draft  layers,  each  with  its  own  cross  sectional  size,  depending  on  the  iceberg 
classification.  These  areas  are  affected  by  the  geostrophic  ocean  water  current  and  a 
calculated  depth  and  time  dependent  local  wind  driven  current.  A  separate  model  is  used 
to  estimate  the  local  wind  driven  current  (Mooney,  1978,  Mountain  and  Mooney,  1979). 
Icebergs  whose  size  are  unknown  are  assumed  to  be  medium  icebergs;  those  whose  shape 
are  unknown  are  assumed  to  be  non-tabular  icebergs.  The  model  is  operated  every  12 
hours  using  the  most  recent  wind  data,  and  drifts  all  icebergs  on  plot  within  the  IIP 
operations  area. 

The  IIP  estimates  that  the  initial  position  error  is  5  nm  regardless  of  sighting 
source  and  that  model  drift  error  increases  linearly  in  5  nautical  mile  per  day  increments 
for  each  24  hours  of  additional  model  drift  up  to  a  maximum  radius  of  30  nautical  miles. 
This  maximum  error  of  30  nautical  miles  occurs  after  5  days  of  drift.  There  is  no  increase 
in  the  maximum  30  nautical  mile  error  estimate  regardless  of  how  long  the  iceberg  is 
drifted  within  the  UP  operations  area.  If  an  iceberg  is  resighted,  the  drift  error  calculation 
is  restarted.  Icebergs  south  of  40°N  are  assumed  to  have  a  daily  drift  error  of  10  nm, 
accumulating  over  a  period  of  5  days  to  a  maximum  error  of  55  nm. 
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Iceberg  Drift  Model. 

The  Iceberg  Drift  Model  (Mountain,  1980)  is  comprised  of  a  set  of  four  differential 
equations  that  balance  iceberg  acceleration  with  the  forces  associated  with  air  and  water 
drag,  the  Coriolis  acceleration  and  a  sea  surface  slope  term  which  describes  the  mean 
ocean  currents.  The  equations  of  motion  for  the  drift  of  an  iceberg  in  component  form 

are; 

^  =  (7  (0 

dt 

±  =  y  (2) 

d/ 

^^fV-fV,  ^{\l2p,CD,AJV^  sin  9+j:\l^>.CD^ME,-U)S,V M  (3) 

dt  *  *=i 

^^fU-fUg+[[l  2p„CD^A„W^  cos  (9+  X 1  /  ^^CD„A,  {N,  -  V)S,  ]  /  M  (4) 


where 


x,y  =  east,  north  components  of  position 

U,F  =  east,  north  components  of  iceberg  velocity 

Ug,Fg  =  east,  north  components  of  mean  (geostrophic)  water  velocity 

^  components  of  Ekman  current  in  /th  layer 

E. '  =  Ug+Ug^ ,  east  component  of  total  water  current  in  /th  layer 

=  Vg+Vg ,  north  component  of  total  water  current  in  /th  layer 
S  =  [(£,  -  C/)^  +  (iV,  -  F)^  ,  speed  of  iceberg  relative  to  /th  layer 

Ag  =  iceberg  cross  sectional  area  above  the  water 

A.  =  iceberg  cross  sectional  area  in  the  /th  layer 

/  =  the  Coriolis  parameter 

M  =  the  iceberg  mass 

W  =  wind  speed 

0  =  wind  direction 

CD„,CD„,  =  drag  coefficient  in  air,  water 
Pa  -  A.  =  density  of  air,  water 


Equations  (l)-(4)  are  solved  using  a  fourth  order  Runge-Kutta  method  to  obtain  the  east 
and  north  components  of  the  iceberg  drift. 


The  Ekman  current  components  are  obtained  using  a  local  wind  driven  current 
model  (Mooney,  1978;  Mountain  and  Mooney,  1979)  that  uses  a  96  hour  wind  history. 
The  IIP  implementation  uses  eight  12  hour  periods.  The  model  from  Mountain  and 
Mooney  (1979)  involves  two  equations  that  are  applied  at  each  layer; 


k-\ 


(5) 
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k-l 

where 

F/c  —  magnitude  of  wind  stress  during  the  Ath  interval 
$1.  =  direction  of  wind  stress  during  the  Ath  interval 

Ki^,<Pk  =  constants  affecting  the  magnitude  and  rotation  of  the  current  in  Ath 
interval 

The  local  wind  driven  current  speed  and  direction  are  given  as: 

Current  speed  =  (7) 

Current  direction  =  tan“* (C/ ^ .  /  Fg- )  (8) 

Mountain  and  Mooney  (1979)  present  selected  values  of  Kk,<Pk  for  sixteen  six- 
hour  time  intervals  (covers  a  96-hour  time  history)  and  selected  latitudes.  The  IIP 
implementation  utilizes  base  values  for  45‘’N  and  adjusts  for  other  latitudes  using  a 
gradient  approach.  The  IIP  implementation  uses  eight  twelve-hour  time  intervals  to  cover 
the  96-hour  time  history  of  local  wind.  Mountain  and  Mooney  do  not  specify  how  the 
wind  stress,  F^,  is  to  be  computed.  In  Mooney's  (1978)  development,  the  wind  stress  is 
linearly  dependent  on  the  wind  speed,  F^  ^PoCdW.  In  a  review  of  the  Mooney  model, 
Dick  (1991)  noted  that  many  oceanographers  use  a  quadratic  relationship,  F^.  =p^Cj^W^. 
Despite  the  linear  model  in  Mooney's  formulation,  the  IIP  implementation  does  use  the 
quadratic  form  to  estimate  wind  stress  (lines  49  and  50  in  SUBROUTINE  NEWWIND). 


QUALITATIVE  AND  EMPIRICAL  EVALUATION  OF  THE  MODEL 

Initial  model  tests  in  1980  used  the  tracks  of  2  large  tabular  icebergs,  a  large 
pinnacle  iceberg,  and  a  freely  drifting  satellite-tracked  buoy  to  compare  the  model 
performance  with  actual  iceberg  drift  (Mountain,  1980).  Results  ranged  from 
approximately  a  5  nautical  mile  error  for  a  3  day  drift  to  a  constant  50-80  nautical  mile 
error  in  the  25  day  case.  The  assumed  cause  for  the  error  in  this  test  was  stated  to  be 
inaccurate  wind  and  current  data  inputs  to  the  model. 

In  1985,  drift  model  tests  were  held  in  several  different  parts  of  the  IIP  operations 
area  (Murphy  and  Anderson,  1985).  Four  case  studies  were  performed  using  the  drift 
model.  In  3  of  the  4  cases,  the  drift  of  the  icebergs  as  depicted  by  the  model  using  the 
system  data  had  location  errors  ranging  from  40  nautical  miles  after  a  2.5  day  drift,  30 
nautical  miles  after  a  3.3  day  drift,  and  45  nautical  miles  after  a  4  day  drift.  These  all 
exceed  the  standard  drift  error  assumed  by  the  IIP  of  a  maximum  of  30  nautical  miles  after 
a  drift  of  5  days.  The  4th  case  study  drift  did  remain  well  within  the  standard  drift  error 
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and  did  not  exceed  an  error  of  more  than  approximately  11  nautical  miles  over  a  4  day 
drift.  Better  performance  results  from  the  model  were  realized  when  using  the  observed 
current  and  wind  data  for  all  four  case  studies,  instead  of  the  automatically  provided 
system  data.  (It  is  not  known  whether  current  from  the  local  wind  driven  current  was 
used  at  the  different  iceberg  layers.)  However,  even  with  using  real  time  data,  in  only  one 
case  was  the  predicted  iceberg  drift  position  error  well  within  the  accepted  limits. 
Projecting  a  drift  experiment  such  as  this  for  a  total  period  of  2  weeks,  or  3  times  as  long 
as  these  case  studies,  suggests  that  the  drift  errors  would  continue  to  generate  and  become 
even  larger. 

Using  the  on  scene  wind  and  current  data  resulted  in  estimated  positions  closer  to 
the  actual  iceberg  position  than  using  geostrophic  currents  and  FNMOC  winds.  The 
limited  experiments  suggest  that  the  structure  of  the  model  is  sound  and  that  its  accuracy 
depends  on  the  accuracy  of  the  input  data.  The  experiments  did  not  isolate  wind  or 
current  as  the  primary  causal  factor  in  generating  errors.  However,  the  case  of  one 
iceberg  in  Murphy  and  Anderson  (1985)  showed  significant  improvement  by  using 
observed  currents  with  FNMOC  winds.  It  should  be  noted  that  three  of  the  cases  in  the 
Murphy  and  Anderson  (1985)  experiment  were  located  in  areas  where  significant 
reductions  in  geostrophic  currents  were  instituted  in  1989-90  (Murphy,  Hanson,  and 
Tuxhom,  1990).  It  is  also  important  to  note  that  based  on  the  trends  in  the  various  drift 
results,  if  the  tests  were  continued  beyond  the  4  days,  the  model  error  would  likely  exceed 
the  30  nautical  miles  maximum  drift  error  after  a  period  of  10  days  of  predictions  even 
with  using  the  best  on  scene,  observed  data  available.  In  practice,  it  is  not  known  to  what 
extent  positional  errors  affect  the  overall  accuracy  of  the  predicted  positions. 


ANALYTICAL  EVALUATION  OF  THE  ICEBERG  DETERIORATION  MODEL 

In  analyzing  the  drift  model,  it  is  appropriate  to  separately  analyze  the  local  wind 
driven  (Ekman)  current  model  separately,  and  then  evaluate  the  overall  drift  model. 

Local  Wind  Driven  Current  Model. 

Analytical  Analysis. 

The  east  component  of  the  local  wind  driven  current  is  given  as  follows: 

Ar-l 

The  input  variables  are  the  wind  speed  and  direction.  The  first  order  sensitivity  of  the 
current  speed  with  respect  to  the  inputs  is  given  by  the  partial  derivatives  as  follows: 

JV  T 

^  =  2p„CZ).«isin(0,-«.,)  (10) 
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and 


(11) 


6d^ 

Similar  partial  derivatives  can  be  obtained  for  the  north  component  of  the  local  wind 
driven  current.  With  respect  to  local  wind  speed,  the  change  in  the  east  and  north 
components  varies  with  the  wind  speed  and  greater  increases  (or  overestimates  of  wind 
speed)  will  have  a  greater  effect.  With  respect  to  wind  direction,  the  effect  of  changes 
varies  as  the  cosine  of  the  angular  error.  Unfortunately,  it  is  virtually  impossible  to 
analytically  determine  what  the  impact  will  be  on  the  resulting  current.  The  changes 
represented  by  equations  (10)  and  (1 1)  are  for  the  Ath  time  interval.  The  total  impact  of 
changes  or  data  errors  over  all  time  intervals  is  incorporated  in  the  sum  in  equation  (9)  for 
the  east  component  and  in  equation  (6)  for  the  north  component.  The  overall  effect  on 
current  speed  and  direction  is  finally  reflected  in  equations  (7)  and  (8).  The  derivative  of 
current  speed  is  easily  obtained  and  indicates  that  it  varies  linearly  with  the  wind  speed  and 
proportional  to  the  cosine  of  the  wind  direction  error.  The  derivative  with  respect  to  wind 
direction  can  also  be  obtained,  but  does  not  provide  any  immediate  insight  into  how  the 
current  direction  is  affected  by  input  data  errors.  Clearly,  this  requires  an  empirical 
sensitivity  analysis  to  identify  these  effects. 

Empirical  Analysis. 

There  are  three  models  for  local  wind  driven  current  that  are  available.  Mooney's 
(1978)  model  uses  a  linear  relationship  for  wind  stress.  It  includes  a  table  of  eight  values 
for  (Ar*.,^jt)  corresponding  to  six-hour  increments  of  wind  information  and  an  eddy 
viscosity  of  50  cmVsec  for  latitudes  ranging  fi-om  0°  to  65°N  in  5°  increments.  Mountain 
and  Mooney's  (1979)  model  does  not  specify  how  wind  stress  is  computed.  It  includes  a 
table  of  sixteen  values  for  corresponding  to  six-hour  increments  of  wind 

information  and  an  eddy  viscosity  of  100  and  500  cmVsec  for  latitudes  at  25°  and  50°N 
for  surface  currents  and  the  current  at  20m.  The  third  model  is  the  one  that  is  coded  in  the 
IIP  computer  software  and  used  in  the  DMPS  model.  The  current  code  is  extracted  fi-om 
SUBROUTINE  NEWWIND  in  Figure  1 . 

The  IIP  drift  model  as  extracted  from  Figure  1  and  using  the  surface  current 
Ekman  parameters  from  SUBROUTINE  NEWWIND  at  45°N  latitude  was  modeled  in  an 
Excel  spreadsheet.  Nominal  wind  speed  was  20  knots  and  the  wind  was  from  the  south 
(180°).  Assuming  that  this  wind  was  constant  over  96  hours,  the  resulting  local  wind 
driven  current  at  the  surface  is  10.287  cm/sec  (0.2  knots)  at  086°  using  the  IIP  drift 
model.  The  detailed  results  are  included  in  Appendix  I. 

The  primary  focus  here  is  the  effect  on  the  local  wind  driven  current  when  there 
are  errors  or  perturbations  in  the  input  wind  speed  and  direction.  These  effects  were 
examined  jointly  considering  a  10%  change  in  wind  speed  (±  2  knots)  and  a  15°  change  in 
wind  direction.  The  effect  on  current  speed  is  indicated  in  Table  1. 
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DO  20  1=1,4 
EKX(I)  =  0.0 
EKY(I)  =  0.0 
20  CONTINUE  . 

IN  =  ITIME  -  7 
DO  40  JT  =  IN, ITIME 
ILAG  =  JT  -  IN  +  1 

WNDSPD=IWDATA(IW,JW,  JT)/1 00 

WINDDIR=(IWDATA(IW,JW,JT)-IWDATA(IW,JW,JT)/1 00)*1 00)*1 0 
WNDSQUAR=WNDSPD*WNDSPD*2642.0 

STRESS=ADC*AIRDNSTY*WINDSQUAR  (ADC=0015;  AIRDNSTY’^.OOIZS) 
DO  30L=1,ISIZE 

ROTA  =  ROT(ILAG,L)+(ALAT-45.0)*GRDR(ILAG,L) 

DEK  =  DK(ILAG,L)+(ALAT-45,0)*GRDK(ILAG,L) 

WIND  =  (WINDDIR+ROTA)*0.01 7452925 
EKX(L)  =  EKX(L)  +  DEK*STRESS*DSIN(WIND) 

EKY(L)  =  EKY(L)  +  DEK*STRESS*DCOS(WIND) 

30  CONTINUE 

40  CONTINUE _ 


Figure  1.  Local  wind  driven  current  coding,  SUBROUTINE  NEWWIND. 


Table  1.  Effect  on  current  speed  with  wind  speed  and  direction  perturbations. 


Current  Speed  (cm/sec) 

Wind  direction  (degrees) 

Wind  speed  (kts) 

165 

180 

195 

%  Change 

18 

8.333 

8.333 

8.333 

-19% 

20 

10.287 

10.287 

10.287 

- 

22 

12.477 

12.447 

12.447 

21% 

As  expected  from  the  analytical  results,  changes  in  wind  direction  do  not  have  any 
effect  on  current  speed.  Changes  in  wind  speed,  however,  do  have  an  effect  on  current 
speed.  In  this  case,  a  10%  change  in  wind  speed  has  approximately  a  20%  change  in 
current  speed.  The  same  results  hold  if  the  nominal  wind  speed  is  increased  to  40  knots 
and  a  10%  perturbation  is  effected. 

A  percentage  perturbation  makes  little  sense  with  respect  to  direction.  Instead,  the 
effects  of  a  ±  15°  change  was  examined.  The  results  are  included  in  Table  2. 


As  indicated  in  the  analytical  section,  there  is  little  insight  provided  by  the 
equations.  The  empirical  results  demonstrate  that  changes  in  the  wind  speed  do  not  affect 
the  current  direction  when  the  wind  is  constant  over  the  period.  (The  W  term  factors  out 
of  the  east  and  north  components  and  does  not  affect  the  arctangent  computation.)  What 
is  interesting  in  this  analysis  is  how  the  current  direction  changes.  In  particular  changing 
from  a  wind  direction  of  180°  to  195°  results  in  a  165°  shift  in  the  current  direction.  The 
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15°  decrease  in  the  wind  direction  results  in  a  15°  decrease  in  current  direction.  At  this 
point,  there  is  no  obvious  explanation  for  this  counterintuitive  result. 


Table  2.  Effect  on  current  speed  with  wind  speed  and  direction  perturbations. 


Current  Direction  (degrees) 


Wind  direction  (degrees) 


Wind  speed  (kts) 

165 

180 

195 

18 

71.2 

86.2 

-78.8 

20 

71.2 

86.2 

-78.8 

22 

71.2 

86.2 

-78.8 

%  change 

-17% 

- 

-191% 

Dick's  (1991)  Evaluation. 

Recently,  Dick  (1991)  prepared  an  Interim  Report  on  an  ongoing  evaluation  of  the 
Mooney  (1978)  model  as  used  for  Search  and  Rescue  (SAR)  planning.  She  compared  the 
SAR-Mooney  model  with  a  one-dimensional  version  of  the  mixed  layer  model  developed 
by  Mellor  and  Yamada  (1982)  and  found  that  the  SAR-Mooney  model  was  inconsistent. 
For  comparison  purposes,  the  IIP  drift  model  was  exercised  and  the  comparative  results 
are  in  Table  3. 


Table  3.  Comparative  drift  model  results. 


Mellor-Yamada 

IIP-Mooney  | 

Wind  Speed 
m/s 

Current 

m/s 

Direction 

Current 

m/s 

Direction 

Current 

m/s 

Direction 

2.5 

0.01 

76 

0.04 

48 

0.006 

86 

5.0 

0.03 

76 

0.08 

48 

0.022 

86 

10.0 

0.08 

59 

0.15 

48 

0.097 

86 

20.0 

0.20 

49 

0.30 

48 

0.389 

86 

40.0 

0.47 

43 

0.60 

48 

1.557 

86 

The  results  in  Table  3  indicate  that  none  of  the  models  agree.  For  low  wind 
speeds,  the  IIP -Mooney  model  is  much  closer  to  the  Mellor-Yamada  model,  but  at  high 
wind  speeds,  it  diverges  rapidly.  A  very  interesting  result  is  the  difference  between  the 
SAR-Mooney  model  and  the  IIP-Mooney  current  model.  Apparently,  there  are  two 
different  implementations. 

Dick(1991)  included  some  preliminary  recommendations  (pending  completion  of 
the  analysis)  that  the  SAR-Mooney  model  be  replaced  by  a  more  accurate  model  for  SAR 
and  iceberg  drift  purposes.  If,  in  fact,  the  Mellor-Yamada  model  is  state-of-the-art,  the 
above  comparison  suggests  that  the  IIP-Mooney  model  may  be  overestimating  local  wind 
driven  current  velocity.  This  will  increase  the  inaccuracy  of  the  iceberg  drift,  but  it  is 
impossible  to  identify  the  direction  of  the  error. 
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Summary. 


The  local  wind  driven  current  model  is  an  important  element  of  the  iceberg  drift 
model.  The  analytical  and  empirical  analyses  indicate  that  the  current  speed  magnitude  is 
somewhat  sensitive  to  errors  in  the  wind  speed  that  is  provided  as  an  input.  Similar 
analyses  with  respect  to  errors  in  wind  direction  indicate  that  errors  in  wind  direction  may 
have  a  significant  impact  on  estimated  current  direction.  Other  comparisons  with  a  SAR 
current  model  suggest  that  it  should  be  replaced.  The  IIP  current  model,  based  on  the 
same  foundation  as  the  SAR  current  model,  yields  significantly  different  results.  Before 
any  current  model  changes  are  made  based  on  the  analysis  of  the  SAR  current  model,  the 
differences  between  the  two  models  must  be  resolved. 

Main  Drift  Model. 

The  main  variables/parameters  that  affect  the  drift  model  in  equations  (l)-(4)  are 
environmental  parameters  (wind  speed  and  direction),  geostrophic  current,  and  iceberg 
size  and  shape.  The  interaction  among  these  parameters  determines  how  the  estimated 
drift  will  vary. 

Environmental  Parameters. 

The  sensitivity  results  from  the  local  wind  driven  current  model  obviously  affect 
the  solution  of  equations  (l)-(4).  The  wind  speed  and  direction  also  impact  equations  (3) 
and  (4)  with  respect  to  the  above  water  area  of  the  iceberg.  The  partial  derivatives  of  this 
portion  of  the  equation  will  be  identical  to  equations  (lO)-(ll)  with  the  water  parameters 
replaced  by  the  corresponding  air  parameters.  In  equations  (3)  and  (4),  the  drag 
coefficient  of  the  iceberg  in  air  is  0.0015  and  in  water  is  0.8  (in  the  latest  version  of  the 
code).  The  density  of  air  is  0.00125.  Assuming  that  the  wind  speed  is  ten  times  as  great 
as  the  current  speed,  the  weighting  factor  for  the  air  wind  stress  is  (.00125 *.00 15*  100  = 
.0001875)  compared  with  water  current  stress  value  of  (1*0.8*  1  =  0.8).  The  impact  of 
the  water  force  is  over  4,000  times  as  great  as  the  air  stress,  assuming  equal  areas.  Based 
on  this  analysis,  the  environmental  parameters,  except  for  the  development  of  the  local 
wind  driven  current,  have  little  direct  impact  on  iceberg  drift. 

Iceberg  Classification. 

The  iceberg  classification  regarding  size  (growler,  small,  medium,  large)  and  type 
(tabular,  pinnacled)  determines  the  surface  areas  and  the  iceberg  mass  used  in  equations 
(3)  and  (4).  Mountain  (1980)  provided  values  for  these  parameters  (see  Table  1  in 
Mountain  (1980)).  Those  values  are  currently  used  with  the  exception  of  growlers. 
Current  parameters  for  growlers  is  M  =  9.6x10^ kg,  =  57m^,  and  =  414m^. 

To  explore  the  impact  of  iceberg  size  selection  on  drift  model  results,  we  have 
computed  the  ratio  of  the  iceberg  area  (above  the  water  and  the  four  segments  below  the 
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water)  to  the  iceberg  mass  (4  /  M).  This  factor  appears  in  each  term  in  equations  (3)  and 
(4).  As  discussed  above,  assume  that  the  forces  on  the  dry  area  have  negligible  impact 
because  of  the  density  differences.  Further,  assume  that  the  current  at  all  levels  is  the 
same  (speed  and  direction)  and  is  equal  to  one  speed  unit.  Then  the  "current  contribution" 
in  equations  (3)  and  (4)  is  given  in  the  last  column  of  Table  4. 

Table.  4.  Iceberg  Area/Mass  ratio. 


Dry 

0-20m 

20-50m 

50-1 00m 

100- 

120m 

Current  Contribution 

Non-Tabular 

Growler 

5.94 

43.13 

0.00 

0.00 

0.00 

43.13 

Small 

3.07 

10.40 

10.93 

0.00 

0.00 

21.33 

Medium 

1.01 

2.00 

2.11 

3.00 

0.00 

7.11 

Large 

0.36 

0.64 

0.68 

0.96 

0.25 

2.54 

Tabular 

Small 

2.65 

7.76 

10.61 

0.00 

0.00 

18.37 

Medium 

1.24 

2.03 

2.72 

4.01 

0.00 

8.76 

Large 

0.63 

0.87 

1.18 

1.75 

0.61 

4.41 

The  "current  contributions"  in  Table  4  contribute  to  the  acceleration  of  the  iceberg 
and  when  equations  (l)-(4)  are  solved  will  result  in  those  with  greater  accelerations  having 
greater  velocity.  The  values  in  Table  4  satisfy  our  intuition  that  other  things  being  equal, 
smaller  icebergs  will  drift  faster  than  larger  icebergs.  The  data  in  Table  4  seem  to  suggest 
that  larger  tabular  icebergs  will  drift  faster  than  corresponding  non-tabular  icebergs 
whereas  smaller  non-tabular  icebergs  will  drift  faster  than  smaller  tabular  icebergs.  The 
present  policy  assumes  that  an  iceberg  is  a  medium,  non-tabular  iceberg  if  a  positive 
classification  can  not  be  made.  Although  this  appears  to  be  a  conservative  policy,  the 
"current  contribution"  is  almost  three  times  as  great  for  the  medium  iceberg  as  compared 
with  the  large  (non-tabular)  iceberg.  None  of  the  experiments  to  date  have  attempted  to 
identify  possible  effects  of  misclassification. 

Geostrophic  Currents  and  Positional  Accuracy. 

Finally,  the  last  input  that  impacts  the  solution  of  equations  (l)-(4)  is  the 
geostrophic  current.  Because  the  geostrophic  current  is  an  average  of  past  observations, 
it  is  inherently  accurate  as  to  its  intended  representation.  The  degree  to  which  it 
corresponds  to  actual  current  is  unknown  in  real  time  except  for  those  cases  where  drift 
buoys  are  available.  In  fact,  the  drift  buoys  are  used  to  provide  a  replacement  for  the 
geostrophic  currents  when  real  time  drift  buoy  current  data  is  available.  Accuracy  in 
geostrophic  data  used  in  the  drift  model  is  very  dependent  on  positional  accuracy. 

To  examine  the  impact  of  positional  accuracy,  we  would  like  to  estimate  the 
probability  that  an  iceberg  is  actually  located  in  the  area  for  which  a  geostrophic  current  is 
selected.  The  geostrophic  current  file  is  developed  on  a  20  second  grid.  Assume  that  the 
geostrophic  current  in  adjacent  north/south  grids  is  approximately  the  same,  but  that 
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east/west  grids  may  have  significant  differences,  particularly  with  regard  to  current  speed. 
The  nP  assumes  that  the  initial  error  in  sighting  an  iceberg  is  10  nm  on  the  first  day, 
increasing  by  5  nm  per  day  up  to  a  maximum  error  of  30  nm.  This  error  distribution  is 
normally  represented  as  a  bivariate  normal  distribution  and  the  maximum  error 
corresponds  to  3o  the  range  fi'om  -3g  to  +3g  covers  99.7%  of  possible  locations.  Under 
these  conditions,  we  can  assume  that  the  marginal  density  of  location  across  lines  of 
longitude  is  normally  distributed,  and  with  the  maximum  error  of  30  nm,  c  =  10  mn.  In 
the  following,  we  develop  a  general  model  that  can  be  used  to  estimate  the  probability  of 
interest. 

Assume  that  the  grid  size  is  d  and  that  a  perpendicular  across  the  grid  is  divided 
into  n  segments,  each  of  length  d/n  as  illustrated  in  Figure  2.  Assume  further,  that  an 
iceberg  position  can  be  selected  at  any  point  in  a  given  segment  i.  For  ease  of 
computation,  we  assume  that  this  position  is  at  the  mid-point  of  the  segment.  Its  value  is 
easily  computed  as  fi  =  {2i-\)d  l  {2n)  where  i  is  the  segment  number  (1  to  n).  Assume 
that  the  error  in  position  is  normally  distributed  with  mean  equal  to  the  mid-point  of  an 
arbitrary  segment  and  standard  deviation  equal  to  o. 


Figure  2.  Estimation  of  positional  accuracy. 

The  endpoints  of  the  interval  are  easily  computed  as  Zj^=-{2i-\)dHlna)  and 
Za={d-(S2i-\)dl{2n))y a  and  the  conditional  probability  is  easily  obtained  from 
tables..  The  resulting  probability  is  conditional  with  respect  to  the  given  segment. 
Assuming  that  each  segment  in  the  interval  of  width  d  is  equally  likely,  the  unconditional 
probability  that  the  iceberg  is  in  the  given  interval  of  width  d  is  the  simple  average  of  the 
segment  probabilities. 

In  the  present  case,  r/  =  20  nm  and  a  =30  nm.  With  ±30  nm  error,  an  easy  estimate 
of  the  probability  is  1/3.  Applying  the  above  model,  yields  a  probability  of  0.61  that  given 
a  position  of  an  iceberg  is  estimated  in  a  given  interval,  it  is  actually  there  .  This  amounts 
to  the  probability  that  the  correct  geostrophic  current  is  selected.  The  results  for  various 
allowed  drift  errors  are  included  in  Table  5. 

Table  5  provides  some  insight  into  the  potential  benefits  of  improving  the  position 
estimation  of  the  icebergs.  In  particular,  if  the  initial  position  is  much  more  accurate,  for 
example,  within  3  nm,  the  probability  that  one  will  select  the  correct  geostrophic  current 
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increases  to  0.97.  With  the  correct  current,  there  is  a  much  higher  probability  that  the 
position  at  the  next  update  will  be  correct  and  that  the  correct  values  of  the  geostrophic 
current  will  be  used  in  equations  (3)  and  (4). 

Table  5.  Probability  of  selecting  correct  geostrophic  current  using  temporal  error 

estimates. 


Day 

Error 

Probability 

1 

10  nm 

.87 

2 

15  nm 

.80 

3 

20  nm 

.74 

4 

25  nm 

.67 

5 

30  nm 

.61 

Summary. 

The  analytical  evaluation  of  the  iceberg  drift  model  reveals  that  there  is  little  need 
for  improved  estimates  of  the  environmental  parameters  for  direct  use  in  the  drift  model. 
The  analysis  illustrated  the  importance  of  correct  classification.  It  is  suggested  that  the 
current  policy  of  classifying  unknown  icebergs  as  non-tabular  medium  icebergs  be 
reexamined.  Finally,  the  sensitivity  of  positional  accuracy  was  clearly  illustrated  for  the 
drift  model.  It  is  important  to  be  able  to  improve  the  initial  positioning  accuracy  to  ensure 
that  the  probability  of  using  the  correct  geostrophic  current  is  maximized. 


SUMMARY  AND  CONCLUSIONS 

This  evaluation  and  analysis  of  the  IIP  iceberg  drift  model  has  concluded  that  the 
model  appears  to  be  a  very  reasonable  representation  of  the  drift  process.  This  conclusion 
is  supported  by  our  review  and  limited  empirical  data.  However,  it  appears  that  the 
submodel  providing  the  local  wind  driven  current  should  be  revisited.  The  inconsistencies 
in  the  drift  results  using  three  supposedly  comparable  models  strongly  suggests  that  the 
reasons  for  the  differences  be  explained.  The  policy  of  characterizing  unidentified  icebergs 
as  non-tabular  medium  icebergs  should  also  be  revisited.  The  total  area  to  mass  ratio 
suggests  that  a  slower  drift  may  be  desirable.  Such  a  review  has  to  be  taken  in 
conjunction  with  other  policy  implications  for  the  model  (e.g.,  deterioration  model).  The 
analysis  highlighted  the  critical  importance  of  improving  the  accuracy  of  estimated  iceberg 
positions.  The  flow  through  of  this  source  of  uncertainty  should  be  examined  in  a 
simulation  model  when  the  local  wind  driven  current  model  questions  are  resolved. 
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Appendix  L  Local  Wind  Driven  Current  Data  Analysis. 

This  Appendix  includes  an  analysis  of  the  sensitivity  of  the  local  wind  driven  current  speed 
and  direction  in  the  face  of  perturbations  of  the  values  of  wind  speed  and  wind  direction. 
The  models  implement  the  current  version  of  the  local  wind  driven  current  model 
developed  by  Mooney  (1978)  and  implemented  by  the  International  Ice  Patrol. 
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Surface  Layer 
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Figure  I-l.  Local  wind  driven  current,  10%  speed  perturbation  at  180°direction. 
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Figure  1-2.  Local  wind  driven  current,  10%  speed  perturbation  at  -15°direction. 
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Figure  1-3.  Local  wind  driven  current,  10%  speed  perturbation  at  +15°direction. 
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(Ikl  =  0.5144444  m/s-above  result  for  wind  speed  of  10  m/s) 


Figure  1-4.  Local  wind  driven  current;  comparison  with  Dick  (1991). 
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